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Numerical  tests  of  possible  models  for  the  oxygen  dependence  of  the  divacancy  introduction  rate 
in  silicon  electron  irradiated  at  room  temperature  were  performed  on  a  computer.  Only  the  model 
in  which  oxygen  traps  Si  self-interstitials  can  reproduce  all  the  experimental  data.  Our  modeling 
results  (in  conjunction  with  the  experimental  data)  imply  that  during  room-temperature  electron 
irradiation  of  Si  the  indirect  production  of  di  vacancies  can  be  more  important  than  the  direct 
production  of  V via  the  single-collision  process. 

PACS  numbers:  61.70.Yq,  61.80.  —  x,  61.70.  —  r,  66.30.Lw 
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I.  INTRODUCTION 

In  a  previous  paper  (Part  I)  we  reported  experimental 
results  about  the  influence  of  the  interstitial  oxygen  content 
of  a  Czochralski  (Cz)  silicon  wafer  on  the  divacancy  produc¬ 
tion  during  room-temperature  electron  irradiation.1  We 
found  that  the  divacancy  [V2]  introduction  rate  varied  lin¬ 
early  with  interstitial  oxygen  concentration.  Similarly  the 
introduction  rate  of  vacancy-oxygen  centers  (VO)  displayed 
a  linear  dependence  on  interstitial  oxygen  content.  The  most 
remarkable  result  of  our  experiment  was  that  the  di  vacancy 
concentration  correlated  much  better  to  the  density  of  intro¬ 
duced  vacancy-oxygen  centers  than  to  the  initially  present 
dissolved  oxygen  concentration  and  that  the  relationship 
between  V2  and  VO  was  perfectly  linear 

The  VO  density  of  a  Si  sample  after  electron  irradiation 
is  an  indicator  of  how  many  single  vacancies  were  available 
for  capture  by  the  oxygen  atoms,  i.e.,  the  VO  concentration 
is  proportional  to  the  steady-state  concentration  of  vacan¬ 
cies  and  the  length  of  the  irradiation  (neglecting  saturation 
and  other  effects).  The  linear  relationship  between  V2  and 
VO  indicates  a  possible  dependence  of  the  V2  concentration 
on  the  steady-state  vacancy  concentration,  i.e.,  the  impor¬ 
tance  of  the  indirect  production  of  divacancies. 

Based  on  our  results  we  proposed  different  possible 
models  for  the  increase  in  divacancy  introduction  rate  with 
increasing  oxygen  content.  It  is  the  purpose  of  the  present 
paper  to  test  these  models  numerically.  Of  the  possible  diva¬ 
cancy  enhancement  mechanisms  mentioned  in  Part  I  only 
one  can  reproduce  all  our  experimental  results. 

II.  DEFECT  PRODUCTION  MODELING 
A.  Basle  model 

In  our  computer  model  of  defect  production2"4  we  as¬ 
sume  that  single  vacancies  V  and  Si  self-interstitials  I  are 
created  at  a  constant  rate  G0  during  room-temperature  elec¬ 
tron  bombardment  of  silicon.  Simultaneously  the  irradiation 
introduces  immobile  divacancies  V2  and  a  corresponding 
number  of  interstitials  at  a  constant  rate  GVl0  •  The  mobile 

"'Permanent  address:  Industrial  Research  Institute  for  Electronics  HIKI, 

H-I392,  Budapest.  Hungary. 

"'Permanent  address:  Forsvarets  Forskingsanstalt,  S-58III,  Linkoping, 

Sweden. 


interstitials  and  vacancies  diffuse  and  can  annihilate  each 
other  or  they  can  be  captured  by  trapping  centers.  In  addi¬ 
tion  two  vacancies  are  allowed  to  combine  and  form  a  diva¬ 
cancy.  The  capture  rate  between  any  two  defects  is  assumed 
to  be  of  the  form  4zr(Z>l  +  D2)RI2,  where  Z>,  and  D2  are  the 
diffusion  coefficients  of  the  defects  and  Rl2  is  an  effective 
capture  radius.5-6  For  the  boron-doped  Cz  silicon  crystal 
ip  s  712  cm)  which  we  used  in  our  experiment1  we  restrict 
ourselves  to  these  impurities:  Varying  interstitial  oxygen 
concentrations,  a  substitutional  carbon  density  of  5.0  X  1016 
at/cm\  and  a  substitutional  boron  content  of  2.0X  1015  at/ 
cm3.1 

It  is  known  that  the  interstitial  oxygen  0,  captures  va¬ 
cancies  to  form  vacancy-oxygen  centers  VO,7  while  the  sub¬ 
stitutional  carbon  C,  and  boron  B,  impurities  trap  Si  self¬ 
interstitials.8-9  Interstitial  carbon  C,  is  mobile  at  room 
temperature  and  can  react  with  other  defects.  Interstitial  bo¬ 
ron  B,  complexes  with  other  Si  self-interstitials  to  form  large 
interstitial  clusters 

In  the  simplest  version  of  our  model  the  following  de¬ 
fect  reactions  are  assumed: 


V  +  I—*d>  (recombination), 

(1) 

v+  V^V2, 

(2) 

Vi  +  I-*V, 

(3) 

V  +  O/— >vo, 

(4) 

vo  +  /-+o„ 

(5) 

I  +  cs^c„ 

(6) 

c,  +  V-+Cs, 

(7) 

VO  +  C,— *0,0,, 

(8) 

B,  +  I— »B,, 

(9) 

B,  +  /— »B,/, 

(10) 

B,/  +  /^B,/2, 

(ID 

B,/„  +  /— *B,/„ ,  , . 
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TABLE  I.  Simultaneous  differential  equations  for  the  reactions  (1-11). 


S[V\ 


‘G0  —  R  '(V,  IW,  +  Dy)[V][I\  -  1R  '(K  V\Dr[VV 

m 

-  R  ‘(V,  0,\Dy[V)V>,]  -  R '1C,.  VI DCl  +  Dy){C, ][KJ 

+  R{V„I\D,[vm 

=Go  +  2Gyi0-R'[V,IW,+Dy)[VUI] 


- 11  '(VO, /iO,(VO][/]  -  R  ‘[Vv  7tf>,[VJ(7] 
-  R  '(C„  7)A[C*J[/]  -  R  '(B,,7)(B,](7] 


a  iv,} 


ii  =  6.,o  +  2R  'fV,  V]Dy[V\2  -  R  '(K2, 7)7),  [  V,]  (7 ) 


a[0 


=  -  ff  '(O,.  V\Dy{Q,\[V\  +  R  '(VO,  7)7>,[VOM7] 
at  % 

=  R  '(O,.  V]Dy{0,\[V]  -  R  (VO,  7)7),[VO][7] 
-ff'(C„V010C([C(][V01 


=  - 11 '(CJt,7>DJ[Cs]l7]+ 11  (C„KHZ)c, +/)„)[€,]  [  KJ 
at 

=  R  '(Cj,  I]D,  [Cs  ]  (7 1  -  R  '(C„  KH7>c,  +  ^(C,  ]}V] 
-ff'(C„VO)7>C([C,][VO] 

=  —  H  '(Bj,7)7),[Bs](7] 

c»[B,l 


dt 


=  R  ’(Bs,  7|7>/(BS][7] 


always  11 '  =  (X,  Y)  =  4 irR  (JT,  Y |. 


These  reactions  require  the  system  of  simultaneous  differen¬ 
tial  equations  shown  in  Table  I."  We  solved  this  system  nu¬ 
merically  by  using  a  computer.  Values  for  the  various  phys¬ 
ical  parameters  which  we  used  in  the  calculation  are  shown 
in  Table  II. 

All  capture  radii  are  taken  to  be  5  A  which  is  appropri¬ 
ate  for  neutral  defects  and  elastic  interaction  at  300  K. 12  We 
neglect  charge  state  effects,  although  they  could  be  impor¬ 
tant.12  The  values  of  the  diffusion  coefficients  for  the  differ¬ 
ent  species  which  are  shown  in  Table  II  are  based  mainly  on 


TABLE  II.  Input  parameters  used  in  numerical  computations. 


Capture  radii:  R{X,  Y )  =  5  A  for  all  defects  X  and  T. 

Diffusion  coefficients  ;300  K); 

D,  =  3.16x10  *cm7s. 

Dy  =4.I5X  10  "cm Vs, 

Dc,  =  1.0x10  14  cm  Vs, 

=  1.0x10  i:  cm;/s. 

Generation  rates. 

G„  =5.0xl0'6  cm  '/s, 

Gyp  =  2.0x  10'*cmVs. 

Initial  values  of  defect  concentrations: 

(O,]  =  varying:  5.0x  10'"atotns/cm'  to  1. Ox  10"  atoms/cm', 
[C*}  =  5.0x  10“  atoms/cm', 

[8*1  =  2.0x  I015  atoms/cm’. 

All  others  =  0. 


experimental  values  for  the  migration  energy  and  the  as¬ 
sumption  that  the  jump  frequency  is  equal  to  the  Debye  fre¬ 
quency.13  Our  confidence  in  D„  is  good.13*17  However,  we 
are  much  less  certain  about  the  accuracy  of  D,,'6'7  DC/,'* 
and  Ionization  effects  on  the  diffusion  of  the  defects 
are  neglected  in  this  paper.  The  vacancy-interstitial  genera¬ 
tion  rate  shown  in  Table  II  is  converted  into  fluence  rate  by 
assuming  that  the  number  of  V,  I  pairs  produced  per  2-MeV 
electron  is  0.520  and  that  only  10%  of  all  V,  1  pairs  escape 
correlated  recovery.  Table  II  also  shows  the  initial  concen¬ 
trations  of  the  various  defects. 


B.  Application  to  the  question  of  the  role  of  oxygen  in 
di  vacancy  production 

In  Part  I  we  proposed  different  possible  mechanisms  for 
the  oxygen  dependence  of  the  di  vacancy  production  rate  in 
electron-irradiated  silicon  at  room  temperature.  We  tested 
those  suggestions  numerically  with  our  computer  model  of 
radiation  damage.  Below  we  list  the  different  mechanisms 
mentioned  in  Part  I  and  the  necessary  modifications  to  the 
model  described  in  Sec.  II  A. 

1.  Kinetics 

Only  the  solution  of  the  system  shown  in  Table  I  is 
required. 

2.  Breiot’s  Suggestion 21 

We  have  to  include  the  capture  of  Si  self-interstitials  by 
interstitial  oxygen  and  the  possible  interactions  of  the  result¬ 
ing  complex  with  other  defects.  The  following  defect  reac¬ 
tions  were  added  to  Eqs.  (1)— (1 1): 

O,  +  I— >0,1,  (12) 

0,1  +  V—*0,.  (13) 

3.  VgO  Center 

We  allow  here  for  the  formation  of  a  V20  complex  from 
VO.  The  respective  concentrations  of  V20  and  V2  are  added 
to  get  an  "apparent”  divacancy  concentration.  We  added  the 
following  reactions  to  Eqs.  ( 1 H 1 1 ): 

V  +  VO— vV20,  (14) 

V20  +  /— VO.  (15) 

4.  Impurity  dependence  of  primary  damage 

a.  Generation  rate  ofV,  I  pairs.  The  introduction  rate  of 
V,  I  pairs  is  assumed  to  depend  on  the  interstitial  oxygen 
content.  We  therefore  write 

G(0;)  =  Go+C[O,].  (16) 

b.  Direct  generation  rate  of  divacancies.  We  allow  the 
direct  generation  rate  of  V2  to  depend  on  the  interstitial  oxy¬ 
gen  concentration.  The  constant  GVr  is  therefore  replaced  by 

G„,(0,)  =  GVj0+C[O,].  (17) 

For  both  cases  (a  and  b )  we  used  C  =  0. 1  s_l. 
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Interstitial  Oxygen  Concentration  UO^cm'3) 


FIG.  1.  Numerical  results  for  the  divacancy  concentrations  vs  increasing 
initial  oxygen  concentration  after  a  total  electron  dose  of  !.0x  10'"  elec- 
t rons/cm'.  "Brelot,"  “V.O  +  V.,“  “G  =  G|0,|,”  ”G,.  =  GK|0,|,"  and 
"Kinetics”  denote  the  models  described  in  Ref  (1  j  and  in  Sec.  II  B 

III.  RESULTS 

In  Fig.  1  we  show  the  divacancy  concentrations  which 
we  obtain  after  a  dose  of  l.Ox  1018  electrons/cm2  for  the 
different  models  described  in  Sec.  II  B  versus  varying  initial 
interstitial  oxygen  contents.  Only  Brelot’s  model  and 
Gy  (O, )  give  increasing  [  V2]  versus  increasing  [O, .]  All  other 
models  give  either  nearly  constant  [F2]  (kinetics,  i.e.,  the 
simplest  system)  or  decreasing  [K2]  versus  increasing  inter¬ 
stitial  oxygen  [V20  +  V2  or  G(0,)]. 

Figure  2  shows  the  vacancy-oxygen  concentrations  ob¬ 
tained  in  the  same  calculations.  All  five  models  show  an  in- 


Interstitial  Oxygen  Concentration  (lO^etn  3) 


FIG.  2.  Numerical  results  for  the  vacancy-oxygen  concentrations  vs  in¬ 
creasing  initial  oxygen  concentration  obtained  in  the  same  calculations. 


Voconcy-Oiygen  Concentration  ( I0IT cm*3) 


FIG.  3.  The  numerical  data  of  Figs.  I  and  2,  but  here  plotting  divacancy 
concentration  vs  the  corresponding  vacancy-oxygen  concentration. 

creasing  A-center  density  versus  rising  initial  amounts  of 
oxygen.  In  Brelot’s  model  there  exists  a  nearly  perfect  linear 
relationship  between  [VO]  and  [  O,  ]  over  the  whole  range  of 
[O,  ]  values.  This  is  the  case  for  the  other  models  only  for 
larger  [O,]. 

The  relationship  between  [  V2]  and  [VO]  after  an  elec¬ 
tron  dose  of  1. Ox  1018  electrons/cm2  is  shown  in  Fig.  3.  The 
models  V..O  +  V2,G  (O, )  and  the  simplest  system  (kinetics) 
all  show  declining  divacancy  concentrations  with  larger  va¬ 
cancy-oxygen  concentrations.  In  Brelot’s  model  [F;]  rises 
nearly  linearly  with  [VO],  while  for  Gv  (O,)  the  divacancy 
concentration  is  also  increasing,  but  the  relationship  to  [VO] 
is  not  linear. 

These  qualitative  trends  of  the  different  models  are  sta¬ 
ble  with  respect  to  variation  of  the  input  parameters,  i.e., 
diffusion  coefficients,  capture  radii,  and  C  (in  Eqs.  16  and 
17),  within  reasonable  limits.  We  therefore  believe  that  the 
results  displayed  in  Figs.  1-3  reflect  intrinsic  characteristics 
of  the  models  described  in  Part  I  and  Secs.  II  B  1— II  B  4,  and 
they  are  not  due  to  the  fortuitous  choice  of  a  particular  set  of 
parameters. 

We  did  not  attempt  to  fit  the  experimental  data  quanti¬ 
tatively,  since  (as  pointed  out  in  Part  I)  it  was  not  possible  to 
convert  the  measured  absorption  coefficient  values  into  de¬ 
fect  concentrations. 

IV.  DISCUSSION 

Our  experimental  findings  reported  in  Part  I  give  us 
three  criteria  which  have  to  be  met  by  the  curves  shown  in 
Figs.  1-3.  As  mentioned  in  Sec.  I  we  found  that  both  [  F2]  and 
[VO]  increased  linearly  with  increasing  interstitial  oxygen 
concentration  and  [F2]  vs  [VO]  was  also  linear.  A  look  at 
Figs.  1-3  shows  that  only  the  numerical  results  of  Brelot’s 
model  compare  favorably  with  these  results.  Only  this  model 
reproduces  qualitatively  the  three  linear  relationships  found 
by  our  experiments.  The  other  models  are  in  qualitative  dis¬ 
agreement  with  our  experimental  data.  On  the  basis  of  this 
we  believe  that  the  formation  of  0,1  pairs  is  the  cause  for  the 
increase  in  the  divacancy  introduction  rate  observed  form- 
creasing  interstitial  oxygen  concentrations. 
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This  conclusion  is  consistent  with  the  annealing  behav¬ 
ior  of  the  0,1  pair,  since  the  available  experimental  evidence 
indicates  that  the  0,1  pair  has  the  thermal  stability  required 
for  a  room-temperature  irradiation.22  Although  the  0,1  -as¬ 
sociated  IR  band  at  935  cm  ~ 1  is  largest  after  low-tempera¬ 
ture  (~77  K)  irradiation  and  decreases  upon  annealing  to 
310  K  it  is  still  present  even  after  annealing  at  350  K.22  We 
therefore  conclude  that  0,-1  interaction  can  take  place  dur¬ 
ing  room  temperature  electron  irradiation  and  that  the  re¬ 
sulting  complex  can  be  stable  at  room  temperature. 

Further  support  for  the  correctness  of  our  conclusion 
comes  from  the  following  consideration.  In  Brelot’s  model 
the  divacancy  increase  in  oxygen  containing  silicon  is  not 
related  to  the  chemical  nature  of  the  oxygen  impurity,  but 
only  to  the  fact  that  oxygen  functions  as  a  Si  self-interstitial 
trap.  Therefore  replacing  oxygen  by  any  other  interstitial 
trap  should  effect  a  similar  divacancy  increase.  This  actually 
has  been  observed  for  carbon.23,24  Alternatively  we  could 
directly  vary  the  steady-state  vacancy  concentration  by 
changing  the  concentration  of  vacancy  sinks  (impurities, 
surfaces,  etc.),  and  we  should  observe  a  divacancy  change. 
Such  an  experiment  has  been  done  by  Wang  et  al.2S  They 
studied  the  surface-defect  distributions  of  electron-irradiat¬ 
ed  silicon  by  transient  capacitance  techniques.  They  ob¬ 
served  a  depletion  of  defects  near  the  surface  and  attributed 
it  to  vacancy  trapping  at  the  surface,  further  they  found  that 
both  VO  and  V2  had  the  same  spatial  dependence.  We  reana¬ 
lyzed  their  data  for  the  divacancy  and  the  VO  center.  Plot¬ 
ting  [FJ  vs  [VO]  gives  the  graph  shown  in  Fig.  4.  We  note 
that  [Fj  depends  linearly  on  [VO].  The  behavior  is  in  quali¬ 
tative  agreement  with  the  curve  of  Brelot’s  model  in  Fig.  3. 
This  further  supports  the  notion  that  oxygen  can  influence 
the  steady-state  vacancy  concentration  by  interstitial  trap¬ 
ping. 

The  significance  of  our  experimental  result  that  [V2] 
can  be  proportional  to  [VO]  (i.e.,  that  [  V2]  reflects  the  steady- 


Voconcy  -  Oxygon  Concentration  (cm'5) 


FIG.  4.  The  date  of  Wang  era/.,”  but  here  plotting  divacancy  concentration 
tE,  -  0.23  eV  energy  level)  vs  vacancy-oxygen  concentration  (£,  -0.18eV 
energy  level). 
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state  single-vacancy  concentration)  and  the  subsequent  ex¬ 
planation  of  this  relationship  by  Brelot’s  model  is  in  its  rela¬ 
tion  to  the  question  of  divacancy  production  per  se.  Since 
there  is  conclusive  experimental  evidence  that  divacancies 
can  form  by  the  direct  collision  process  during  electron  irra¬ 
diation  of  silicon26  27  for  a  long  time  the  formation  of  diva¬ 
cancies  via  combination  of  two  vacancies  was  thought  to  be 
negligible.  We  can  now  answer  the  question  of  what  fraction 
of  a  certain  divacancy  concentration  introduced  by  room- 
temperature  electron  irradiation  was  produced  by  a  direct 
collision  process  and  what  fraction  is  due  to  the  agglomera¬ 
tion  of  single  vacancies.  In  our  experiment  we  found1: 

aK,  =  0.57aVo  +  0-25  cm- '.  (18) 

Since  the  fraction  of  a  K  which  is  proportional  to  avo  is  due 
to  the  combination  of  two  vacancies  while  the  constant  gives 
an  upper  limit  of  directly  introduced  divacancies  we  can 
rewrite  Eq.  (18)  as 

&  V*  ^indirect  4*  ^direct  (19) 

with  ain<Ureet  =0.57aVo  and  aainct  =0.25  cm-1.  For  the 
largest  ayo  we  had  aVi  =  1.2  cm-1.1  We  then  get 
^indirect  :adirec>  =  0.95:0.25  and  we  therefore  have  that  the  in¬ 
direct  production  rate  of  divacancies  can  be  four  times  as 
large  as  the  direct  production  rate  during  room-temperature 
irradiation  with  2-MeV  electrons.  It  is  interesting  that 
Barnes  found  a  similar  ratio  in  a  low-temperature  (76  K) 
neutron-damage  experiment.  Barnes  observed  immediately 
after  irradiation  at  76  K  (i.e.,  with  the  single  vacancies  im¬ 
mobile)  only  25%  of  the  maximum  [V2]  value  which  was 
obtained  after  a  330-K  anneal.28  Therefore  in  his  case  the 
indirect  production  rate  of  V2  was  three  times  as  large  as  the 
direct  production  rate. 

V.  SUMMARY 

Numerical  tests  of  possible  models  for  the  oxygen  de¬ 
pendence  of  the  di  vacancy  introduction  rate  in  at  room-tem¬ 
perature  electron-irradiated  silicon  were  performed  on  a 
computer.  Only  one  of  the  models  can  reproduce  all  the  ex¬ 
perimental  results.'  The  other  models  give  results  which  are 
in  qualitative  disagreement  with  the  experimental  data.  We 
envision  the  mechanism  responsible  for  the  enhancement  of 
divacancy  production  as  follows:  Interstitial  oxygen  pre¬ 
vents  by  the  capture  of  Si  self-interstitials  interstitial-va¬ 
cancy  recombination.  This  in  turn  increases  the  steady-state 
vacancy  concentration  and  ultimately  enhances  the  diva¬ 
cancy  concentration  via  agglomeration  of  two  single  vacan¬ 
cies. 

Our  experimental  and  modeling  results  show  that  the 
indirect  formation  of  di  vacancies  during  room-temperature 
irradiation  of  silicon  with  2-MeV  electrons  can  be  more  im¬ 
portant  than  the  production  of  V2  by  the  single  collision  pro¬ 
cess. 
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